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Abstract 

If the activity of active galactic nuclei (AGNs) is predominantly induced by major 
galaxy mergers, then a significant fraction of AGNs should harbor binary massive 
black holes in their centers. We study the mass function of binary massive black holes 
in nearby AGNs based on the observed AGN black-hole mass function and theory of 
evolution of binary massive black holes interacting with a massive circumbinary disk 
in the framework of coevolution of massive black holes and their host galaxies. The 
circumbinary disk is assumed to be steady, axisymmetric, geometrically thin, self- 
regulated, self-gravitating but non- fragmenting with a fraction of Eddington accretion 
rate, which is typically one tenth of Eddington value. The timescale of orbital decay 
is then estimated as ~ 10 8 yr for equal mass black-hole, being independent of the black 
hole mass, semi-major axis, and viscosity parameter but dependent on the black-hole 
mass ratio, Eddington ratio, and mass-to-energy conversion efficiency. This makes 
it possible for any binary massive black holes to merge within a Hubble time by 
the binary-disk interaction. We find that (1.8 ±0.6%) for the equal mass ratio and 
(1.6 ±0.4%) for the one-tenth mass ratio of the total number of nearby AGNs have 
close binary massive black holes with orbital period less than ten years in their centers, 
detectable with on-going highly sensitive X-ray monitors such as Monitor of All-sky 
X-ray Image and/or Swift/Buist Alert Telescope. Assuming that all binary massive 
black holes have the equal mass ratio, about 20% of AGNs with black hole masses of 
10 6 ' 5_7 M Q has the close binaries and thus provides the best chance to detect them. 

Key words: black hole physics - accretion, accretion disks - binaries:general - 
galaxies:nuclei 
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1. Introduction 



Most galaxies are thought to have massive black holes at their centers (Kormendy & 
Richstone 1995). Massive black holes play an important role not only in the activities of 
active galactic nuclei (AGNs) and quasars but also in the formation and evolution of galax- 
ies(Magorrian et al. 1998; Ferrarese&Merritt 2000; Gebhardt et al. 2000). Galaxy merger leads 
to the mass inflow to the central region by tidal interactions and then a nucleus of the merged 
galaxy is activated and black hole grows by gas accretion (Yu & Tremaine 2002). At some step, 
the outflow from the central black hole sweeps away the surrounding gas and quenches the star 
formation and black hole growth. This also produces an observed correlation between the black 
hole mass and velocity dispersion of individual galaxies (Di Matteo et al. 2005). 

During a sequence of processes, binary massive black holes with a subparsec-scale sepa- 
ration are inevitably formed before two black holes merge by emitting gravitational radiation. 
Recent hydrodynamic simulations showed the rapid binary black hole formation in the parsec 
scale within several Gyrs by the interaction between the black holes and the surrounding stars 
and gas in gas- rich galaxy merger (Dotti et al. 2007; Mayer et al. 2007). Even if there are 
transiently triple massive black holes in a galactic nucleus, the system finally settles down to 
the formation of binary massive black holes by merging of two black hole or by ejecting one 
black hole from the system via a gravitational slingshot (Iwasawa et al. 2006). 

In coalescing process of two massive black hole, there has been the so-called final parsec 
problem: it is still unknown how binary massive black holes evolve after its semi-major axis 
reached to the subparsec scale where the dynamical friction with the neighboring stars is no 
longer effective. Many authors have tackled the final parsec problem in the context of the 
interaction between the black holes and the stars, but there has been still extensive discussions 
(Begelman et al. 1980; Makino 1997; Quinlan & Hernquist 1997; Milosavljevic & Merritt 2003; 
Sesana et al. 2007; Matsubayashi et al. 2007; Matsui & Habe 2009). There is other possible 
way to extract the energy and the angular momentum from binary massive black holes by the 
interaction between the black holes and the gas surrounding them. This kind of the binary- 
disk interaction could also be the candidate to resolve the final parsec problem (Ivanov et al. 
1999; Gould & Rix 2000; Armitage & Natarajan 2002; Armitage & Natarajan 2005; Escala et 
al. 2005; Hayasaki 2009; Cuadra et al. 2009; Haiman et al. 2009) in spite of a claim(Lodato et 
al. 2009). Some authors showed that there exist close binary massive black holes with a short 
orbital period less than ten years and significant orbital eccentricity(Armitage & Natarajan 
2005; Hayasaki 2009; Cuadra et al. 2009). 

Hayasaki et al. (2007) studied the accretion flow from a circumbinary disk onto binary 
massive black holes, using a smoothed particle hydrodynamics (SPH) code. They found that 
mass transfer occurs from the circumbinary disk to each black hole. The mass accretion rate 
significantly depends on the binary orbital phase in eccentric binaries, whereas it shows little 
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variation with orbital phase in circular binaries. Periodic behaviors of the mass accretion rate 
in the binary system with the different geometries or system parameters were also discussed 
by some other authors(Bogdanovic et al. 2008; MacFadyen & Milosavljevic 2008; Cuadra et 
al. 2009). Recently, Bogdanovic et al. (2009) and Dotti et al. (2009) proposed the hypothesis 
that SDSSJ0927f 2.65+294344.0 consists of two massive black holes in binary, by interpreting 
the observed emission line features as those arising from the mass-transfer stream from the 
circumbinary disk. 

Hayasaki et al. (2008) have, furthermore, performed a new set of simulations at higher 
resolution with an energy equation based on the blackbody assumption, adopting the same 
set of binary orbital parameters we had previously used (Hayasaki et al. 2007) (a = 0.01 pc, 
eccentricity e = 0.5, and mass ratio q = 1.0). By this two-stage simulation, they found that 
while the Optical/NIR light curve exhibits little variation, the X-ray/UV light curve shows 
significant orbital modulation in the triple-disk system, which consists of an accretion disk 
around each black hole and a circumbinary disk around them. X-ray/UV periodic light variation 
are originated from a phase-dependent mass transfer from circumbinary disk (cf. Hayasaki 
&; Okazaki 2005). The one-armed spiral wave on the accretion disk induced by the phase- 
dependent mass transfer causes the mass to accrete onto each black hole within one orbital 
period. This is repeated every binary orbit. These unique light curves are, therefore, expected 
to be one of observational signatures of binary massive black holes. 

Highly sensitive X-ray monitors over a wide area provide us with a unique opportunity 
to discover close binary massive black holes in an unbiased manner, based on the detection of 
the orbital flux modulation. Monitor of All-sky X-ray Image (MAXI; Matsuoka et al. 2009), a 
Japanese experimental module attached to the International Space Station, is now successful 
in operation since the launch in 2009 July. MAXI, covering the energy band of 0.5-30 keV, 
achieves a significantly improved sensitivity as an all X-ray monitor compared with previous 
missions. According to the hard X-ray luminosity function of AGNs by Ueda et al. (2003), 
~1,300 nearby AGNs can be detected at the confusion flux limit of ~0.2 mCrab from the 
extragalactic sky at galactic latitudes higher than 10°. Among them, the brightest ~100 AGNs 
can be monitored every 2 months with a flux accuracy of 20% level. Over the plausible mission 
life of MAXI (> 2yr), it is possible to detect binary massive black holes in nearby AGNs with 
binary orbital periods less than 10 yr. Besides MAXI, the Swift/Burst Alert Telescope (BAT) 
survey (Tueller et al. 2010) can make a similar job in the hard X-ray band of 15-200 keV. 

In this paper, we investigate mass functions of binary massive black holes with a very- 
short orbital period detectable with MAXI and/or Swift/BAT. The plan of this paper is orga- 
nized as follows. In Section 2, we describe the evolutionary scenario of binary massive black 
holes in the framework of coevolution of massive black holes and their host galaxies. Section 3 
shows mass functions of close binary massive black holes. They can be written as the product of 
the observed black-hole mass function of nearby AGNs and probability for finding binary mas- 
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sive black holes, based on the evolutionary scenario as described in Section 2. Brief discussions 
and conclusions are summarized in Section 4. 

2. Final-parsec evolution of binary massive black holes 

We first describe the evolution of binary massive black holes, focusing on interaction 
with surrounding gaseous disks in the framework of coevolution of massive black holes and 
their host galaxies. 

Does the black hole mass correlate with the velocity dispersion of bulge in individual 
galaxies despite that there is a single or binary in their center? This is one of fundamental 
problems in the framework of the coevolution of massive black holes and their host galaxies. 
In some elliptical galaxies, there is a core with the outer steep brightness and inner shallow 
brightness. Binary massive black holes are considered to be closely associated with such a core 
structure with the mass of stellar light deficit (Ebisuzaki et al. 1991; Milosavljevic & Merritt 
2001). Recently, Kormendy & Bender (2009) showed the tight correlations among black hole 
masses, velocity dispersions of host galaxies, and masses of stellar light deficits, using the 
observational data of 11 elliptical galaxies with cores. This suggests that these correlations are 
still held for even if there are not only a single massive black hole but also binary massive black 
holes in cores of elliptical galaxies. 

Binary massive black holes are considered mainly to evolve via three stages (Begelman 
et al. 1980). Firstly, each of black holes sinks independently towards the center of the common 
gravitational potential due to the dynamical friction with neighboring stars. If the binary can 
be regarded as a single black hole, its gravitational influence radius to the field stars is defined 
as 



where the tight correlation between the black hole mass and one-dimensional velocity dispersion, 
(T* , of the stars, the so-called M bh — a* relation: M bh /1O 7 M = @i (cr*/200kms~ 1 ) /?2 is made use 
of. Unless otherwise noted, (3\ = 16.6 and $2 = 4.86 are adopted by (Merritt & Milosavljevic 
2005) in what follows. 

When the separation between two black holes becomes less than 1 pc or so, angular 
momentum loss by the dynamical friction slows down due to the loss-cone effect and a massive 
hard binary is formed. This is the second stage. The binary harden at the radius where the 
kinetic energy per unit mass of the star with a*, equals to the binding energy per unit mass of 
the binary(Quinlan 1996). Its hardening radius is defined as 




(1) 





where q is the black-hole mass ratio. 
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Finally, the semi-major axis of the binary decreases the radius at which the gravitational 
radiation dominates, and then a pair of black holes merge into a single supermassive black 
hole. The detailed timescale in each evolutionary phase will be described in the following three 
subsections. 

2.1. Star driven phase 

2.1.1. The dynamical friction 

Each black hole sinks into the common center of mass due to the dynamical friction with 
ambient field stars. The merger rate of two black holes is given by(Binney & Tremaine 1987) 
ait) 0.428 , A GM bh 1 

a{t)~ V2 ** a?{ty {6) 

where a(t) is the separation between two black holes and In A « 10 is the Coulomb logarithm. 
The decaying timescale of black-hole orbits is then written 



tdi 



a(t) 



a J \10 7 M Q/ 

where = lOOpc is the typical core radius of the host galaxy. The integration of merger rate 
gives the following equation 

a(t) 




where 



a 



tf ~4.2xl0 6 [yr](^M ( 1 . (6) 

[y 1 v l0( W v 10 M ©/ 

Recall that 02 ~ 5, and hence the index of mass dependence is —4/5 for both t^i and t^ { . 
2.1.2. Stellar scattering 

Even after hardening of the binary, the binary orbit continues to decay by loss cone 
refilling of stars due to the two-body relaxation. In addition, the repeated gravitational slingshot 
interactions with stars makes the orbital decay significantly rapid for black hole with mass less 
than fewxl0 6 M Q . For the system with a singular isothermal sphere, the timescale is given 
as(Milosavljevic & Merritt 2003) 



(5) 



1//32-1 



a(t) 



a(t) 



3.0 x 10V] ( ~T7T | { • (7) 



a(t) J \ W 7 M Q/ 

For the black hole with mass greater than 1O 6,5 M , this mechanism contributes inefficiently 
to the orbital decay of the binary on the subparsec scale. Instead, the binary-disk interaction 
is likely to be a dominant mechanism of the orbital decay. Quite recent N-body simulations 
show that stellar dynamics alone can also resolve the final parsec problem (e.g., Berentzen et 
al. 2009). 
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2.2. Gaseous-disk driven phase 

The circumbinary disk would be formed inside the gravitational influence radius after 
hardening of the binary. The inner edge of circumbinary disk is then defined as 

1\ 2/3 




1-2/02 



(8) 



where a(t) is the semi-major axis of binary, and m = 2 and I = 1 are adopted unless otherwise 
noted (Artymowicz & Lubow 1994). 

For simplicity, the circumbinary disk is assumed to be a steady, axisymmetric, and 
geometrically thin with a differential rotation and fraction of Eddington accretion rate: 



TjM] 



Edd 



2.2 x 10" 



Me 
yr 




(9) 



where rj, e and M^dd are the Eddington ratio, mass-to-energy conversion efficiency, and 
Eddington accretion rate defined by MEdd = (l/e)47rGM bh m p /ccTT, where m p , c, and or show 
the proton mass, light velocity, and Thomson scattering cross section, respectively. 
The surface density of circumbinary disk can be then written as 



2ttu 



dlnr 



(10) 



where the eddy viscosity is defined as v = ac cv H with the Shakura-Sunyaev viscosity parame- 
ter(Shakura & Sunyaev 1973), a, characteristic velocity, c cv , and disk scale-hight, H . 

The stability criterion for self-gravitation of the disk is defined by the Toomre Q-value: 
Q = Vtc CY /iiGYi. If the disk structure obeys a standard diskwitYi M acc with 77 = 0.1, Q is much less 
than 1. This means that the disk is massive enough to be gravitationally unstable. Therefore, 
we introduce a self-regulated, self-gravitating disk model (Mineshige & Umemura 1996; Bertin 
1997). The condition of self-regulated disk is given by 

Q~l. (11) 

From equations (10) and (11), the effective sound velocity of self-gravitating disk is written as 

1/3 

(12) 



GM acc 


lnr 


2a sg 


lnfi 



where we adopt a = a sg ^0.06 and c cv = c sg . If the radiative cooling in the disk is so efficient, the 
disk would fragment. The criterion whether the disk fragment or not is given by a sg = 0.06(Rice 
et al. 2005). If a sg > 0.06, the fragmentation occurs in the disk and causes the subsequent star 
formation. Such a situation is beyond the scope of our disk model. Assuming that the disk 
face radiates as a blackbody, the sound velocity can be written 
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(13) 



y fi I \ aixr A o J 

where R g , fi, and a are the gas constant, molecular weight, and Stefan-Bolzmann constant, 
respectively. 

The self-gravity of the disk is stronger than the gravity of the central black hole at the 
self-gravitating radius where the total disk mass equals to the black hole mass. The rotation 
velocity of the disk then become flat outside the self-gravitating radius: 

-2/3 

- 2/3 m- 2/ v^y /3 (14) 

where we put c cv = c sg . Inside r sg , the angular frequency of the disk corresponds to Keplerian 
one where the gravity of central black hole dominates the dynamics of the disk. When c cv = c s , 
the disk transits from the self-regulated, self-gravitating disk to standard disk. Its radius is 
given as 

N 4/3 /GM acc \~ 8/9 /3GM bh M acc \ 1/3 




stsg 



2a sg J \ &7T(j 
a s „ \ 8 / 9 /0.1 \ / n \ I Muu \ _2/9 



~ 4.4 x 10- 4 [pc] ^ — — MSH • (15) 

As the binary evolves, the disk structure gets to depend on black hole mass. Since n n (ah) is 
less than r sg and more than r stsg in the all black-hole mass range, the circumbinary disk is 
initially modeled as the self-regulated, self-gravitating disk with the Keplerian rotation. When 
the semi-major axis decays to the decoupling radius defined by equations (25) and (26), r in (ad) 
is less than r stsg for 1O 5 M < M bh < 3 x 10 6 M o . The disk structure for 10 5 < M bh < 3 x 1O 6 M 
can then be described by the standard disk theory, whereas the disk still remains to be the 
self-regulated, self-gravitating disk with Keplerian rotation for other black-hole mass ranges 
(see the dotted line of Fig. 2). 

The circumbinary disk and binary exchanges the energy and angular momentum through 
the tidal/resonant interaction. For moderate orbital-eccentricity range, the torque of binary 
potential dominantly acts on the circumbinary disk at the 1:3 outer Lindblad resonance radius 
where the binary torque is balanced with the viscous torque (Artymowicz & Lubow 1994). 

On the other hand, the circumbinary disk deforms to be elliptical by the tidal/resonant 
interaction. The density of gas is locally enhanced by the gravitational potential at the closest 
two points from each black hole on the inner edge of circumbinary disk (Hayasaki et al. 2007). 
The angular momentum of gas is removal by the locally enhanced viscosity, and thus the gas 
overflows from the two points to the central binary. An accretion disk is then formed around 
each black hole by the transferred gas(Hayasaki et al. 2008). The mass transfer therefore adds 
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its angular momentum to the binary via two accretion disks(Hayasaki 2009). 

In a steady state, the mass transfer rate equals to the accretion rate, M acc . Since it is 
much smaller than the critical transfer rate defined by equation (41) of Hayasaki (2009), the 
effect of torque by the mass transfer torque can be neglected. The orbital-decay rate is then 
approximately written by equation (40) of Hayasaki (2009) as 

where e is the orbital eccentricity of the binary and the net torque, J c bd, from the binary to 
circumbinary disk can be approximately written as 

^^(il^-L^* (17) 

q t vis - m M bh VI - e 2 

where t v is,in = r \J v- m = \J GM\^r- m / («c 2 v ) is the viscous timescale measured at the inner edge of 
disk and M d is the local disk mass defined as M ld = 7rrf n £ in . From equation (10), the product 
of the viscous timescale and ratio of the black hole mass to local disk mass is given by 



M ld \ _1 M bh d\nVt 



t . . 1U = 9 



M hh J M a 



dlnr 



(18) 



Substituting equation (17) with equation (18) into equation (16), the orbital-decay rate can be 
expressed by 

a(t) 1 



a(t) tf s ' 

where tf as is the characteristic timescale of orbital decay due to the binary-disk interaction: 

t gas z 



(19) 



2 


q M bh 




3 4 / 3 


(l + g) 2 M acc 


dlnr 




3.1xl0 8 [yr]-Al — ) (_ \ (20 ) 



(1 + q) 2 \ T] 

where we adopt Q = Qk- Note that tf as is independent of the black-hole mass, semi-major 
axis, and viscosity parameter, but dependent on the black-hole mass ratio, Eddington ratio, 
and mass-to-energy conversion efficiency. These arise from the assumptions that the disk is ax- 
isymmetric with a fraction of Eddington accretion rate and its angular momentum is outwardly 
transferred by the viscosity of Shakura-Sunyaev type. 
Integrating equation (19), we obtain 

f—hrO- (21) 

Following equations (29) of Hayasaki (2009), the orbital eccentricity increases with time 
in the present disk model. The orbital eccentricity is, however, expected to saturate during 
the disk-driven phase, because the angular momentum of the binary is mainly transferred to 
the circumbinary disk when the binary is at the apastron. The saturation value of orbital 
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eccentricity becomes e = 0.57. This value is estimated by equating the angular frequency at 
the inner edge of circumbinary disk with the orbital frequency at the apastron. 

2.3. Gravitational-wave driven phase 

The merging rate by the emission of gravitational wave can be written by(Peters 1964) 

a{t) 256G 3 M^ h q /(e) 



as 



a(t) 5c 5 a 4 (l + g) 2 (l-e 2 ) 7 / 2 ' 

where /(e) = 1 + 73e 2 /24 + 37e 4 /96. The coalescent timescale is then given as 

4 . n , „\2 



t 



gw 



a(t) 




fa(t)\ 


a(t) 


32 





(22) 



(23) 



c q 

where r§ = 2GM^/c 2 is the Schwarzschild radius. 

The semi-major axis decays to the transition radius, a t where the emission of gravita- 
tional wave is more efficient than the binary-disk interaction. In other words, t gw becomes 
shorter than tf as inside the transition radius. Comparing equation (20) with equation (23), the 
transition radius is defined by 

1/4 

(24) 



at 
rs 



32 ctf* 



q 



/(e) 



5 r s (l + g) 2 (l-e 2 ) 7 / 2 _ 
When the timescale of orbital decay by the emission of gravitational radiation is shorter 
than the viscous timescale measured at the inner edge, the circumbinary disk is decoupled with 
the binary. The decoupling radius is then defined by 

32 ct viSiS q /(e) 



rs 



5 r s (l + g) 2 (l-e 2 ) 7 /2 
for the standard disk, and 



32 ct 



vis.S 



q 



/(e) 



2/7 



(25) 



(26) 



5 r s (l + g) 2 (l-e 2 ) 7 / 2 _ 
for the self-regulated, self-gravitating disk. Here ty- 1Bj s is the viscous timescale measured at 
the inner edge of circumbinary disk when a(t) = rs- Depending on the disk model, it can be 
written as 

I V V— Y /4 ( Mbh V /4 

1 Wl) V1O 7 M J 

for the standard disk, and 

1/3 /n i \ 2/3 ^ x 2/3 / - , v 1/3 



tvis,s~5.9x 10 2 [yr] 



— ) 



'or 



(27) 



t vis ,s~ 5.6 x 10 4 [yr] 



0.06 



21) f^-Y /3 ( Mbh 

n ) vo.i/ vio 7 m 



a sp \0.1 
for the self-regulated, self-gravitating disk, respectively. 
Integrating equation (22), we obtain 



(28) 
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(29) 





o 



) 



-3 



(1 + g) 2 



(30) 



where eo = 0.57 is the initial orbital eccentricity at a t . 
2.4- Observable period range for binary black holes 

The resonant /tidal interaction causes the mass transfer from the circumbinary disk to 
accretion disk around each black hole(Hayasaki et al. 2007). The ram pressure by mass transfer 
acts on the outer edge of accretion disk and gives a one-armed oscillation on the disk (cf. 
Hayasaki & Okazaki 2005). The one-armed wave propagates from the outer edge to the black 
hole, which allows gas to accrete onto the black hole within the orbital period. This is repeated 
every binary orbit. This mechanism therefore originates periodic light variations synchronized 
with the orbital period(Hayasaki et al. 2008). 

For the observational purpose, a±o is defined as the semi-major axis corresponding to 
feasible orbital period, lOyr, detectable with MAXI and/or Swift/BAT by 



Fig. 1 shows the mass dependence of each orbital period evaluated at a t , a d , and a w for equal- 
mass binary massive black holes. The dashed line, dotted line and horizontal dash-dotted line 
show the orbital period at a t , and aio, respectively. The area filled in with the solid line 
shows the existential region of binary black hole candidates with periodic light-curve signatures 
detectable with MAXI and/or Swift/BAT. 

Fig. 2 shows the mass dependence on each characteristic semi-major axis, a^, a t , a^, and 
aio, for the equal-mass binary. All of the semi- major axes get longer as black hole mass is more 
massive. The decoupling radius is described by equation (25) when Mbh ^3x 1O 6 M , whereas 
it is described by equation (26) when Mbh ^3x 1O 6 M . Note that a 10 is on the track, where 
the binary evolves by the binary-disk interaction, when Mbh^2 x 1O 7 M , whereas aio is on the 
track, where the binary evolves by the emission of gravitational wave, when Mbh^2 x 10 7 M Q . 

Fig. 3 shows the orbital-decay timescale, \a/a\, of binary massive black holes with Mbh = 
1O 7,5 M in panel (a) and corresponding elapsed time in panel (b). The dashed line shows the 
timescale in the first evolutionary phase where the orbit decays by the dynamical friction 
(hereafter, dynamical- friction driven phase). The solid line shows the timescale in the second 
evolutionary phase where the orbit decays by the binary-disk interaction (hereafter, disk-driven 
phase). The dotted line shows the timescale in the final phase where the orbit decays by 
the dissipation due to the emission of the gravitational wave (hereafter, gravitational-wave 




(31) 
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driven phase). The dash-dotted line shows the orbital-decay timescale of by stellar scattering 
(hereafter, stellar-scattering driven phase). We note that the orbital-decay timescale of stellar- 
scattering driven phase is too long for stellar scattering to be efficient mechanism in binary 
evolution. The binary therefore evolves toward coalescence via first dynamical-friction driven 
phase, second disk-driven phase, and final gravitational-wave driven phase. The orbital-decay 
timescale of disk-driven phase is the longest among the other two. 

Fig. 4 shows the orbital-decay timescale of the binary with the same format as that of 
panel (a) of Fig. 3, but for Mbh = 10 6 M Q . Note that stellar scattering is efficient after hardening 
of the binary, as shown in the dash-dotted line of Fig. 4. 
3. Mass function of binary massive black holes 

Observed black-hole mass function of nearby AGNs allows one to study mass functions of 
binary massive black holes based on the evolutionary scenario described in the previous section. 
Fig. 5 shows the fractional mass function of hard X-ray selected AGNs in the local universe, 
where the black hole mass is estimated from the K-band magnitudes of their host galaxies, as 
compiled by Winter et al. (2009). This mass function is 99% complete for the uniform sample of 
local Seyfert galaxies, and hence the uncertainties caused by the sample incompleteness should 
be regarded negligible. 

Assuming that two galaxies start randomly to merge during the interval, £*, where t* is 
the look back time from present universe, a probability for finding binary massive black holes in 
the present universe can be expressed by |o/d|/t*. The number of binary massive black holes, 
iVbbh) in AGNs can be obtained by N hhh = / c (|a/d| /£*)/-/Vagn, where iV&GN is the number of 
AGNs and f c is a dimensionless fraction parameter. More than 25% of the host galaxies of 
Swift/BAT AGNs show evidence for on-going mergers (Koss et al. 2010), where the separation 
between central two cores in the merging galaxies is more than kpc scale. Therefore, we set 
f c = 0.75 unless otherwise noted, since we discuss binary black holes with a separation of much 
smaller scale. 

One can estimate the probability for finding binary black holes in nearby AGNs by 
putting £agn into £*, where £agn = Mbh/Maec shows the e- holding accretion timescale which 
corresponds to the typical lifetime of AGNs. Fig. 6 (a) and (b) show the mass-dependences 
of |a/d|/iAGN, evaluated for q = 1.0 and q = 0.1, respectively. The dashed line, dash-dotted 
line and dotted line show the probability for finding binary black holes with a^, a t , and a^, 
respectively. It is noted from both panels that the probability of is lower than that of a t 
in the mass range less than the order of 1O 6 M , because the binary more rapidly evolves by 
stellar scattering than by disk-binary interaction as shown in the dash-dotted line of Fig. 4. 
The probability evaluated at a t keeps constant over all of mass-ranges. 

The solid line shows the integrated probability for finding binary black holes with the 
semi-major axis less than a w . The integrated probability is approximately given by At/tAGN, 
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Fig. 1. Observable period range of binary massive black holes with characteristic semi- major axes. The 
dashed line shows the mass dependence of the orbital period evaluated at the transition radius, a t , where 
the dominant mechanism of binary evolution changes from the angular momentum loss by the binary-disk 
interaction to the dissipation by emitting gravitational wave radiation. The dotted line shows the one 
evaluated at the decoupling radius, a^, where the circumbinary disk decouples with binary massive black 
holes. The horizontal dash-dotted line shows the orbital period of ten years. The region drawn in solid 
lines shows the mass dependence of the orbital period less than ten years, which is detectable with M AXI 
and/or Swift. 
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At 



tf s ln(a 10 /a t ) + tf 
tr(a 10 )-tr(a d ) 



{a t )-tf(a d ) a 10 >a t 
Oio < at, 



(32) 
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Fig. 2. Mass dependence of characteristic semi-major axes of binary massive black holes. The solid line 
shows the radius corresponding to the orbital period of ten years, Oio- The dashed line and dotted line 
show the transition radius, at, and decoupling radius, ad, respectively. The dash-dotted line shows the 
hardening radius, ah, where the binding energy per unit mass of the binary equals to the kinetic energy 
of a star surrounding the binary. 

from equation (20), (21), (29), and (30). 

The integrated probability estimated for aa < a < a±o is the monotonically decreasing 
function of black hole mass. Note that they rapidly decreases as the black hole mass becomes 
greater than 10 7 M Q for q = 1.0 and 5 x 10 7 M o for q = 0.1. 

Fig. 7 shows mass functions of binary massive black holes in AGNs. The mass function 
is defined by multiplying the black-hole mass function of AGNs by the probability for finding 
binary black holes. The mass functions evaluated at the hardening radius, ah and transition 
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Fig. 3. (a) Orbital-decay timescale, \a/a\, of evolution of binary massive black holes with a semi-major 
axis from lOOpc to 10~ 4 pc. The total black hole mass is M^h = 1O 7 ' 5 M0 with equal mass ratio, q = 1.0. 
(b) Corresponding elapsed time of evolution of binary massive black holes. In both panels, the dashed 
line shows the first evolutionary phase in which two black holes get close each other toward the hardening 
radius, a^, by their angular momentum loss due to the dynamical friction with surrounding field stars. 
The solid line shows the second evolutionary phase from ah to at in which the angular momentum of the 
binary is removed by binary-disk interaction. The binary evolves from to at where the stage changes 
from the disk-driven phase to the gravitational-wave driven phase in which binary massive black holes 
finally coalesce by the emission of the gravitational radiation. There is the decoupling radius, ad, in the 
third evolutionary phase shown in the dotted line. In panel (a), the dash-dotted line shows the timescale 
of orbital decay by stellar scattering. 
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Fig. 4. Same format as panel (a) of Fig. 3, but for Mbh = 10 6 M Q . 

radius, a t , are exhibited in panel (a) and (b), respectively. In both panels, the solid line and 
the dashed line show the mass function with q = 1.0 and that with q = 0.1, respectively. 

There is less population of binaries with mass less than 10 6 5 M Q in panel (a), because 
the binary more rapidly evolves by stellar scattering than by disk-binary interaction in the mass 
range, as shown in the dash-dotted line of Fig. 4. Total fraction of binary black holes over all 
mass ranges are (4.3 ±2%) for q = 0.1 and (13 ±4%) for q = 1.0 in both panels (the quoted 
errors reflect the statistical uncertainties in the Swift/BAT AGN mass function). It is noted 
that from both panels that binary black holes of Mbh = 10 8,5_9 M Q are the most frequent in the 
nearby AGN population. 

Fig. 8 shows the mass functions of binary black holes with a constraint that the orbital 
period is less than ten years in both cases of q = 1.0 and q = 0.1. From the figure, binary black 
holes of M bh = 1O 6 5 " 7 M for q = 1.0 and those of M bh = 1O 7 - 5_8 M for q= 0.1 are the most 
frequent in the nearby AGN population. It is notable that, assuming that all the binaries have 
equal black-hole mass ratio, 18% of AGNs with black hole of 1O 6 ' 5_7 M has binary black holes. 

Total fraction of binary black holes over all mass ranges are (1.8 ±0.6%) for q = 1.0 and 
(1.6 ±0.4%) for q = 0.1. We can therefore observe 15 ~ 27 candidates for 1300 AGNs detectable 
with MAXI, assuming that activities of all nearby AGNs lasts for £agn = Mbh/M acc . Note that 
MAXI covers the softer energy band (2-30 keV) than Swift/BAT (15-200 keV), and hence the 
ratio of type-1 (unabsorbed) AGNs to type-2 (absorbed) AGNs will be higher in the MAXI 
survey (~8:5 based on the model by Ueda et al. 2003) than in the Swift/BAT survey (~1:1, 
Tueller et al. 2010). Here we have referred to the same AGN mass function, however, since 
we do not find statistically significant difference between the observed mass functions of type-1 
and type-2 AGNs based on the Kolmogorov-Smirnov test. 
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4. Summary & Discussion 



We study mass functions of binary massive black holes on the subparsec scale in AGNs 
based on the evolutionary scenario of binary massive black holes with surrounding gaseous disks 
in the framework of coevolution of massive black holes and their host galaxies. 

As a very recent progress in observations of binary massive black holes with the Sloan 
Digital Sky Survey (SDSS), there is a claim that two broad emission line quasars with multiple 
redshift systems are subparsec binary candidates(Borson & Lauer 2009). The temporal vari- 
ations of such the emission lines are attributed to the binary orbital motion(Loeb 2009; Shen 
& Loeb 2009). These can be used as complementary approaches to search for binary massive 
black holes with MAXI and/or Swift/BAT. 

Recently, Volonteri et al. (2009) predicted the fraction of binary quasars at z < 1 based 
on the theoretical scenario for the hierarchical assembly of supermassive black holes in a ACDM 
cosmology. They adopted the merging timescale of binary black holes with a circumbinary disk 
estimated by Haiman et al. (2009), in order to explain the observed paucity of binary quasars 
in the SDSS sample (2 out of 10000; Bogdanovic et al. 2009; Dotti et al. 2009; Borson & Lauer 
2009). For the black hole mass range of ~ 10 8 M Q , which these SDSS quasars likely have, our 
calculation gives a similar merging timescale (~ 10 8 year, independent of mass). Hence, our 
model will also be compatible with the SDSS results when applied to the same cosmological 
model. In a lower mass range, however, we predict a significantly longer merging timescale, by 
a factor of 10 at ~ 1O 7 M , than that by Haiman et al. (2009), which rapidly decreases with the 
decreasing mass. Hence, much larger fractions of subparsec binary black holes are expected in 
our model than in Volonteri et al. (2009) if low mass black-holes are considered. 

Kocsis & Sesana (2010) studied the nHz gravitational wave background generated by 
close binary massive black holes with orbital periods between 0.1-10 years, taking account of 
both the cosmological merger rate and such the binary-disk interaction as the planetary (type 
II) migration(Haiman et al. 2009). The orbital-decay timescale for low black- hole mass binaries 
(A^bh < 1O 7 M ) is much shorter than that of our model. This suggests that little stochastic 
gravitational wave background is attenuated by applying our model for their scenario, because 
the amplitude of gravitational wave background is proportional to the root of ratio of the 
orbital-decay timescale of the disk-driven phase to that of the gravitational-wave driven phase. 

Our main conclusions are summarized as follows: 

(1). Binary massive black holes on the subparsec scale can merge within a Hubble time by 
the interaction with triple disk consisting of an accretion disk around each black hole and 
a circumbinary disk surrounding them. Assuming that the circumbinary disk is steady, 
axisymmetric, geometrically thin, self-regulated, self-gravitating but non-fragmenting with 
a fraction of Eddington accretion rate, its orbital-decay timescale is given by ~ 3.1 x 
10 8 g/(l + g)(0.1/?7)(e/0.1)[yr], where q, rj, and e show the black-hole mass ratio, Eddington 
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ratio, and mass-to-energy conversion efficiency, respectively. 

(2) . Binary black holes of Mbh = 1O 8 ' 5_9 M in the disk-driven phase are the most frequent 

among the AGN population. Assuming that activities of all nearby AGNs lasts for the 
accretion timescale, Mbh/M acc , the total fraction of binaries with the semi-major axis 
evaluated at the hardening radius and transition radius are estimated as (4.3 ± 2%) and 
(13 ± 4%), respectively. 

(3) . Assuming that all binary massive black holes have the equal mass ratio (q = 1.0), ~ 20% 

of AGNs with Mbh = 1O 6,5 ~ 7 M harbor binary black holes with orbital period less than 
ten years in their center. This black-hole mass range therefore provides the best chance 
to find such close binary black holes in AGNs. 

(4) . The total fraction of close binary massive black holes with orbital period less than ten 

years, as is detectable with MAXI and/or Swift/BAT, can be estimated as (1.8 ±0.6%) 
for q = 1.0 and (1.6 ±0.4%) for q = 0.1. 

We thank anonymous referee for the useful comments and suggestions. We also thank 
Mike Koss and Richard Mushotzky for providing us with the latest result on the merging rate 
of Swift/BAT AGNs before publication. KH is grateful to Atsuo T. Okazaki, Takahiro Tanaka, 
and Stanley P. Owocki for helpful discussions. This work has been supported in part by the 
Ministry of Education, Science, Culture, and Sport and Technology (MEXT) through Grant- 
in-Aid for Scientific Research (19740100, 18104003, 21540304, 22340045, 22540243 KH and 
20540230 YU, and 22740120 IN), and by the Grant-in-Aid for the Global COE Program "The 
Next Generation of Physics, Spun from Universality and Emergence" from MEXT of Japan. 

References 

Armitage, P. J., & Natarajan, P. 2002, ApJ, 567, L9 

Armitage, P. J., & Natarajan, P. 2005, ApJ, 634, 921 

Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 

Begelman, M. C, Blandford, R. D., & Rees, M. J. 1980, Nature, 287, 307 

Berentzen, I., Preto, M., Berczik, P., Merritt, D., & Spurzem, R, 2009, ApJ, 695, 455 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press), 428 

Bertin,G. 1997, ApJ, 478, L71 

Bogdanovic, T., Smith, B. D., Sigurdsson, S., & Eracleous, M. 2008, ApJ, 174, 455 
Bogdanovic, T., Eracleous, M., & Sigurdsson, S. 2009, ApJ, 697, 288 
Borson, T. A., & Lauer, T. R. 2009, Nature, 458, 53 

Cuadra, J., Armitage, P. J., Alexander, R. D., & Begelman, M. C. 2009, MNRAS, 393, 1423 
Di Matteo, T., Springel, V., k Hernquist, L. 2005, Nature, 433, 604 
Dotti, M., Colpi, M., Haardt, F., & Mayer, L. 2007, MNRAS, 379, 956 

Dotti, M., Montuori, C, Decarli, R., Volonteri, M., Colpi, M., & Haardt, F. 2009, MNRAS, 398, 73 



17 



30 




c 



] [ 



8 



10 



L0G(MJ [M a ] 



0- 



Fig. 5. Normalized mass function of massive black holes in AGNs, detected with Swift/BAT (15-200 
kcV) given by Winter ct al. (2009). 



18 
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(a) and q=0.1 (b). The probability is defined as the orbital-decay timescale normalized by the accretion 
timescale, |a/d|/tAGN, where tAGN = M^/ M acc . The solid line shows the integrated probability for finding 
binary massive black holes with the semi-major axis from aio to ad- The dashed line, dash-dotted line, 
and dotted line show the probabilities evaluated at ah, at, and ad, respectively. 
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